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Abstract 

The angle and temperature dependence of the photon scattering rate for 
Bose-stimulated atom recoil transitions between occupied states is compared 
to diffraction and incoherent Rayleigh scattering near the Bose-Einstein tran- 
sition for an optically thin trap in the limit of large particle number, N. Each 
of these three processes has a range of angles and temperatures for which it 
dominates over the others by a divergent factor as N — > oo. 
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The quantum statistics of trapped atoms can drastically effect small angle light scattering 
from a cold, degenerate cloud of gas. The scattering of an incident photon is accompanied 
by a transition between trap states of a recoiling atom. At small enough scattering angle 
or momentum transfer, the final atom state may very well already be occupied. For bosons 
there is an enhancement of the transition rate proportional to rif + 1, where rif is the 
number of bosons already in that final trap state. For fermions, transitions to occupied 
states are blocked; the analogous factor is rif — 1. (This Fock basis description is the natural 
one for the discussion of the phenomenon, particularly in the context of gaseous systems.) 
Considerable attention has been given to the consequent impact on resonant line-shape [U 
the scattered photon energy is shifted up and down by the differences in trap-state energies, 
which are then averaged over the contributing processes in the thermal ensemble. However, 
the enhancement of the scattering rate itself and its angular dependence may also prove to 
be dramatic and interesting effects 0,0]. 

Ref. considered Born approximation scattering per unit volume for a plane wave of 
light in a uniform medium of bosonic atoms as a function of angle and temperature. As the 
temperature, T, drops below the Bose-Einstein critical temperature, T c , as a consequence 
of this Bose enhancement the scattering rate increases sharply for angles corresponding to 
momentum transfers d such that the number of particles per {h/S) 3 volume is large. While 
plane waves and infinite, uniform media are physically unattainable, they represent a reason- 
able approximation to a beam that is many wavelengths wide incident on the central region 
of a cloud that is many times yet larger still. This paper addresses another limit: when the 
whole trap is illuminated uniformly by the incident beam. In particular, I investigate the 
angle and T dependence of the scattering rate for a plane wave incident on an isotropic, 
harmonic trap containing N non-interacting atoms. Again, I consider only the Born ap- 
proximation, i.e., a single scattering, which is appropriate to an optically thin trap. (This 
may be realized by a combination of low density and sufficiently small basic cross-section, 
e.g., achieved by sufficient detuning from resonance.) A useful reference level is given by 
the isotropic Rayleigh scattering or fluorescence rate, and comparison with the even brighter 
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coherent diffraction at yet smaller angles is made. In addition to the differential rates per 
solid angle, I also compare the total amounts of scattered light due to the various processes, 
integrated over all angles. 

There are many possible strategies for doing this conceptually quite straightforward 
calculation. I have chosen to focus on the leading, large-iV behavior of each possible con- 
tributing sub-process for temperatures such that T/T c ~ 0(1). In particular, while the 
trap ground state is treated explicitly, the other states are described semi-classically, e.g., 
with a continuous density of states and stationary-phase evaluations of WKB wave function 
overlaps. This allows analytic evaluation (or at least estimates) of the various contributions 
to the total scattering. It proves quite enlightening because the several sub-processes and 
regions of scattering angle scale differently with N at fixed T/T c . An alternative strategy 
of exact numerical evaluation would require considerable detail over a great range in N to 
extract comparable information. 

Laboratory detection of scattered light at very small angles has already been achieved in 
the context of cold atom traps. Dramatic images of condensation into a trap ground state 
have been obtained from coherently diffracted light [|J. Essentially, one diverges all small 
angle light, blocks out the unscattered forward beam, and then focuses the remainder. Both 
diffraction and Bose-enhanced scattering have the properties that they are brighter (per 
unit solid angle) than the isotropic scattering and they do not contribute to heating the gas. 
Hence, they both make possible non-destructive observation by allowing a greatly reduced 
illuminating intensity for a given detected signal strength. However, diffraction is confined 
to smaller and smaller angles with increasing trap size. Thus, for larger traps of the future, 
clean diffractive signals may be harder to extract. In contrast, while dimmer than diffraction 
at the smallest angles, the Bose-enhanced scattering contributes over a much larger range 
of angles — a range that increases with N — and for much of this region is considerably 
brighter than Rayleigh scattering. 

In an infinite, uniform medium, diffraction cannot be addressed. It is of infinite intensity 
but relegated to zero angle. In contrast, the present analysis allows an explicit comparison 
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of diffraction and Bose-stimulated scattering. This calculation also underscores the fact that 
the latter occurs only to the extent that the recoil from a photon can knock an atom from 
one occupied state to another and is proportional to the product of the two occupation 
numbers. As will become apparent, as iV — > oo only a vanishing fraction of the total 
number of particles actively participate in the Bose enhancement. As a consequence, the 
total luminosity from Bose-stimulated scattering, integrated over angle, is significantly less 
than that from diffraction. On the other hand, diffraction reflects the overall particle density, 
and its angular distribution scales inversely as the size of the system. In contrast, the Bose- 
enhanced scattering depends on the momentum of particles in the occupied states. This 
momentum and, consequently, the typical scattering angle grow with increasing temperature. 

This paper is organized as follows. Section I contains the thermal preliminaries. Sec- 
tion II contains the kinematics and the reference Rayleigh scattering. Section III discusses 
diffraction. Bose-stimulated scattering is evaluated in section IV, first for transitions to 
and from the trap ground state and then for transitions between excited states; these are 
the only "new" formulae of this paper, but they make most sense in the context of and in 
comparison with the other material. The qualitative lessons are summarized in section V. 

I. THERMAL PRELIMINARIES 

It is convenient to use the natural units of the three-dimensional isotropic harmonic 
oscillator. In particular, take 

m = ujq = h = 1 , 

where m is the particle mass and ojq is the oscillator fundamental frequency. Also, measure 
temperature T in energy units (i.e., the trap level spacing) so that Boltzmann's constant 
k B = 1. I will use Cartesian coordinates throughout the analysis. So the single particle 
stationary states are labeled by a triplet of integers m = (m x ,m y ,m z ), and the energy of 
such a state can be taken to be e m = m x + m y + m z . (As is evident from the form of e m , m 
is not a vector with respect to spatial rotations.) 
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Divide the total number of particles N a,s N = Nq + N e , where Nq is the expected 
number of particles in the ground state (also written as (no)) and N e is the expected number 
in excited trap states, i.e., the remainder. For iV — > oo, such a system has a critical 
temperature given by 

T c = iV 1 / 3 ^)- 1 / 3 (1) 

where £(3) ~ 1.202 is the Riemann Zeta function |J. (In eq. (|1|), and throughout, T is 
measured in units of the trap level spacing; hence, to see a transition at fixed, finite, physical 
temperature, one must increase the trap size as N — > oo.) For T < T c , N is a macroscopic 
fraction of N: 

AT. = *(l -(I)'). (2) 

This can be deduced from the following evaluation of N e . One begins by noting that the 
expected occupation of any state m is given by 

(rim) 



; ( Em -/i)/T _ i ' 
where fi is the chemical potential. It follows that 

e^ T = 1 + JL . 

Hence, for all excited states we can take /i/T ~ for T < T c . Thus, for T < T c , 

, 00 \m 2 dm 
N, - 1 2 



o 



e m/T _ J 

T 3 C(3) , (3) 



which, by consistency, implies eq. (0). The integral in eq. (0) is obtained by approximating 
the sum over states m, whose degeneracy is (e m + l)(e m + 2)/2 for large N by an integral 
over e m (now called m). One can also deduce that imposing N e — N for T > T c implies 

Vt>t c „, 18C(3)(T-T C ) , o ([T-T^ 2 " 
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The rate for a boson transition from a particular initial state i to a particular final state 
/ is proportional to Tij(n/ + 1), where n, and nf are the actual occupation numbers. It is this 
product which then must be thermally averaged if the system as a whole is characterized by 
an equilibrium temperature. Hence, in principle, one needs to know the two-level correlations 
{riiUf) as well as the occupation expectations, iVj = (rii). For the purposes of the present 
calculations, it suffices to use the replacement (n^n/) ~ NiNf because the corrections to 
this form lead to changes in any final result that are down by 0(1/N). (Here, "final result" 
refers to the result after the thermal average over all possible initial and final iV-atom states 
is performed.) The reasons are discussed at length elsewhere 0, but, briefly, the issues are 
as follows. For i = f ^ 0, the fractional fluctuations in the occupation of any substantially 
occupied state are large (i. e., 0(1)). However, for fixed N, the different z's fluctuate almost 
independently. So, for example, the variance associated with N e vanishes relative to N e as 
N e — > oo. Also, for i ^ /, the fixing of N induces a correlation between occupations, but it is 
likewise negligible as N — ► oo. Finally, for i = f = 0, fixing N ensures that the fluctuations 
in no are microscopic, even for an ideal gas with T < T c . 



II. GENERAL KINEMATICS AND RAYLEIGH SCATTERING 

A photon of momentum kj is incident on a trap. The goal is to estimate the differential 
cross-section da/dfl for scattering the photon to momentum kf directed in the solid angle Q 
(and also the total cross-section, a). In particular, I integrate over the final photon frequency 
(or, equivalently, over kf). In the natural units defined above, the unit of velocity is of order 
the RMS zero-point velocity of an atom in the trap ground state. In these units, the speed of 
light, c, is very large. If Ae is the change in the struck atom's energy in a particular collision 
and, therefore, minus the change in the photon's energy, then kf = ki — Ae/c. Hence, the 
scattering cross-section da/d 3 k.f is approximately proportional to 8{kf — kj). If the final goal 
is angular dependence, integrated over final photon energy, one can use the virtual equality 
of ki and kf to simplify the kinematics. 
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Let the momentum transfer be defined as 8 = k/ — kj. At small angles, using ki — kf, 8 
is proportional to the scattering angle 9 




dQ~Sd8d v /k*, (4) 

where tp is the azimuthal angle. No information is lost by making the trivial simplification of 
ignoring the effects of photon polarization on angular dependence. (At issue is polarization 
in the scattering plane at large angles; if this were of interest, it could easily be accounted for 
correctly.) With this simplification, the photon angular dependence is (^-independent and 
totally determined by the magnitude of the momentum transfer to the atom, 8. Conversion 
to actual laboratory scattering angles involves knowing the incident photon momentum or 
frequency. However, the discussion of atomic trap physics is simplest and clearest in terms 
of 8. 

Born approximation scattering separates naturally into three types of processes: diffrac- 
tion, Bose-stimulated scattering, and Rayleigh scattering. By "diffraction" I mean the col- 
lisions in which the struck atom remains in its original trap state. Hence, the iV-atom final 
state is the same as the iV-atom initial state, irrespective of which atom is struck. There- 
fore, the amplitudes for photon scattering off each of the N atoms must be added coherently, 
and then the sum is squared to get the rate. When the struck atom changes state, even if 
unobserved, the rates for all the different possible processes are added (incoherently). As 
mentioned above, these rates are proportional to rij(n/ + 1). The rii ■ 1 piece is what is 
traditionally identified as Rayleigh scattering. It is independent of the occupation of the 
final state. The • rif piece is what I mean by Bose-stimulated scattering. It is an analog 
for bosonic atoms of stimulated emission for photons. 

In all cases, the key transition matrix element that determines the scattering connecting 
the trap state i to / is (i\e^' r \f) . This can be visualized in momentum space as the overlap 
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of the initial momentum wave function shifted by the momentum transfer 6 with the final 
state momentum wave function. However, many of the calculations are easiest to evaluate 
in position space. 

The problem at hand is separable, e.g., in Cartesian coordinates. Hence, the station- 
ary state wave functions (r\i) are products of three one-dimensional harmonic oscillator 
wave functions. Furthermore, the potential and the thermal distribution are isotropic. For 
convenience, we can choose one of our axes, e. g., x, to point along S, and then 

(m\e i6 - r \m'} = (m x \e iSx \m' x ) 5 my ^ y S mz , m > z . (5) 

The one-dimensional wave functions of interest are the ground state, 

( X \0) = TT-VV* 2 / 2 , 

and the WKB wave functions for m 1. For example, for x > 

. . . /7r\ — V 2 / ox -i/4 f sin ( C \f2m - x 2 dx + f ) for x < a 

(x\m) ~ - (2m -x 2 ) ' { . KJx . rx 4 ( 

V2/ [ i exp (— J a \/x 2 — 2mdx) ioi x > a 

where a = \j2m is the classical turning point. 

Rayleigh scattering serves as a convenient reference comparison for the other processes. 
The Rayleigh differential cross-section is the total Rayleigh cross-section for a single atom 
times the following expression: 
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= N . 



Hence, Rayleigh scattering is independent of T and 6 (i.e., isotropic) and proportional to 
N. 

For the sake of simplifying all formulae and facilitating the comparisons of the various 
processes, I suppress the factor of the one-particle cross-section and write all differential and 
total cross-sections in units of that fundamental area. Hence, I write 

G^Raylcigh _ 
CRayleigh = 4:71 N . 



To reiterate, the displayed "cross-sections" are not areas in the natural trap units; rather, 
they all have a single, common factor, the atomic physics one-particle cross-section, sup- 
pressed. 



III. DIFFRACTION 



The diffractive cross-section is given by 

^""diffraction 



i\e i5 - r \i) 



The i th term in the sum is just the Fourier transform of the position space density of the 



,-th 



one-particle state, weighted by the occupation number. Hence, Born approximation 



diffraction has a very classical interpretation, independent of the particle statistics: it is 
always given by the square of the Fourier transform of the density. (Of course, that density 
and its thermal behavior may themselves be very dependent on quantum statistics.) 

For large iV and T ~ 0(T C ), it is appropriate to single out the ground state and treat 
the excited states semi-classically. In particular, the normalized position-space density for 
classical one-particle states of energy m in a three dimensional isotropic harmonic potential 
is 



7r 2 m 2 



m — 



for r 2 < 2m and otherwise. And the ground state density is just 

Po (r) = tt-s/V 3 . 
So, the large- iV diffractive cross-section is 

d 0"di ffr act ion 



d 3 re iS - r 



N p (r) + 



00 \m 2 dm 



e ( m -^)/ T - 1 



Pm(r) 



o 4 Jo z 6 



A e- 52 /4 + 



4T 
~8* 



for /i — 



(6) 
(7) 
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(The z integral in eq. (|6]) is an approximation to an exact infinite sum that arises in the 
evaluation of the preceding line; the integral representation is valid for 5 2 T ^> 1.) Note that 

4T _ ( 2 V N * 
~¥ ~ C(3) ' 

Also, if the characteristic linear size of the diffuse cloud of atoms in excited trap states at 
temperature T is called L, then L 2 ~ 0(T); the corresponding, characteristic momentum 
transfer supported by diffraction from such a cloud, 5l, would be 0(1/ L); so 8\T ~ 0(1). 
Thus, the two terms in (^diffraction/^ have a clear interpretation. The diffraction cross- 
section is nominally 0(N 2 ); however, the N$ part has support only for 5 < 0(1), i.e., the 
inverse of the spatial extent of the ground state. The excited state part has power fall-off 
in 5 and is only 0(N 2 ) or larger for S < 0(T~ l l 2 ) ~ (^(iV -1 / 6 ), where the latter estimate 
applies for T ~ 0(T C ). Of course, in the square of the amplitudes, there is also a cross term 
between the two. 

[For completeness, I note that for fi < the z integral in eq. @ can be expressed as a 
sum of two "modified Bessel functions of the second kind" of ranks and 1 and argument 
8y/—2fi. However, the qualitative and asymptotic behavior is easiest to extract from the 
integral representation itself.] 

To integrate over angles, one must recall eq. (|j) for dfl. The excited-state part of the 
diffractive cross-section appears to be singular as 5 — > 0. However, that is a reflection of 
the inappropriateness of representing the discrete states as a continuum when considering 
very small shifts in momentum. If one simply excludes the vanishingly small forward cone 
as iV — * oo in which the diffuse cloud diffraction contributes, the remaining diffraction is 
dominantly off the condensate: 

^diffraction — TjN . (8) 

Alternatively, one can estimate where the discreteness of the spectrum would cut off the 
integral. For excited states with energies of 0(T), the difference in RMS momentum from 
one level to the next is 0(T~ 1 ^ 2 ). Using that as an estimate for the smallest meaningful 5 
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in integrating eq. (0), one learns that the integrated cross-section due to excited states alone 
is 0(N^ 3 /kf), which does not compete with eq. (§) for T < T c as N — ► oo. 

The explicit appearance of the incident photon momentum ki is inevitable here (and also 
below) when considering total cross-sections. While the Rayleigh differential cross-section is 
genuinely isotropic, diffraction and Bose-stimulated scattering are only significant at small 
angles. Within those ranges they can be much larger than the Rayleigh rate, as indicated by 
the formulae. However, the extent of those ranges is governed by the ratio of the acceptable 
struck-atom momentum transfer (which depends on trap geometry and temperature) to the 
incident photon momentum, which is typically much larger. 

IV. BOSE-STIMULATED SCATTERING 

The Bose-stimulated rate has two sources: scattering atoms into and out of the conden- 
sate and scattering atoms between excited trap states. For each of these classes of processes, 
one first determines the square of the matrix element of e l ^' r between initial and final states. 
Again, it is convenient to use 8 to define one of the Cartesian coordinate directions, e.g., 
8 = 5x. Then the relevant overlap is an integral over one-dimensional oscillator wave func- 
tions. Also, only states which share the same y and z quantum numbers contribute, as 
indicated in eq. ([5]). So the thermal sum over occupations of possible states involves first a 
projection onto the contributing states and then a sum over the thus-projected occupations 
weighted by the non-trivial one-dimensional matrix-element-squared. 

A. In and out of the condensate 

With the choice x = 8/5, only states labeled by m = (m, 0, 0) have non-zero matrix ele- 
ments of e"^' r with the ground state. The one-dimensional overlap of relevance for transitions 
between the trap ground and excited states is 

F{m,5) = |(0|e ifa |m)| 2 . 
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For large N and T ~ 0(T C ), only m 3> 1 is of interest. F is sharply peaked around 5 V2m, 
the maximum classical momentum of a one-dimensional harmonic oscillator of energy m. 
Its width in 8 is 0(1), which is vanishingly small on the scale of the typical value of \fm. 
Hence, it suffices to represent F as a 5-function, whose normalization can be deduced by 
integrating F over 5. The result is 

5 2 

F(m, 5) ~ 5(m ) . 

2 

This behavior has a simple classical interpretation. To bring an energetic particle in 
a three-dimensional harmonic potential to rest with a single impulse and have it remain 
at rest, the following must be true: the energetic orbit must pass through the origin; the 
impulse must be applied when the particle passes through the origin; and the impulse must 
be precisely that required to bring the particle to rest. 

The Bose-stimulated differential cross-section is obtained by summing over m the equal 
contributions of transitions — > m and m — > 0: 



dQ u 7o e"*/r-r 

2N 



e <5 2 /(2T) _ i 

J 2N e~£ for 5 2 > 2T 
\ AN ^ for 5 2 < 2T 

For 5 ~ (9(1) this is O(N N^ 3 ) while the diffractive cross-section is O(Nq); however, the 

latter then falls off exponentially with 5 2 while the <-> m stimulated rate only falls as 

1/S 2 . Note also that this stimulated rate is much larger than the Rayleigh rate as long as 

5 < 0(T 1 / 2 ). When 5 2 > T, the <-> m stimulated rate falls off exponentially with 5 2 /2T 

because of the fall-off in thermal occupation of appropriately high momentum states. 

This analysis is again inadequate for very small 5, where the discreteness of the spectrum 

is relevant. However, we can integrate down to 5 ~ 0(1) for which these approximations 

are sufficient and estimate 

^>l)-^l°g(2T). 
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B. Between excited states 



The excited states m and m' must be related by m = (m, m y , m z ) and m' = (m', m y , m z ) 
to get a non-zero overlap with e^' r . The relevant one-dimensional matrix element is 



G(m,m',S) = \ {m\e lSx \m!) 



This is to be summed against the thermal occupations projected onto the one dimension: 



V{m,m',T,n) = [(e (m+my+mz ~ M)/T - i)( e K+™«+™*-M)/T _ ^ 

m y ,m z 

= T 2 P 



T 



T 



where 



P(a,b) 



z dz 



, {e z+a - l){e z + b - 1) 
Thus, in terms of the functions so defined, 



(9) 



dn 



dm dm' G(m, m , 5) Vim, m , T, /i) . 



P(a, b) can be expressed in closed form in terms of exponentials, logarithms, and dilog- 
arithms. Again, though, its qualitative and limiting behavior is easiest to see from the 
defining integral, eq. (H). P(a,b) is symmetric in a and b, and it is 0(1) when both a and 
b are themselves 0(1). If b < a <C 1, it is 0(log(l/a)). And it is (9(exp(— a — b)) when 
a, b > 1. (Recall, also, that /i/T is for T < T c and C((T - T c )/T c ) for T > T c .) 

G(m,m' ,5) can be evaluated by stationary phase using the WKB wave functions. For 
each m, m', and 5 (all three positive), there is a single position-space point of stationary phase 
in the overlap integral, corresponding to the unique classical position where an impulse of 
magnitude S can change a particle of energy m to one of energy m! . G(m, m', 5) is symmetric 
in m, m'; the following simple form is for m > m'\ 

G(m,m>,5) = \ 2A 2 \ 

y or if 2m o < [m — m — ^) 



13 



Unlike F(m,5), the support of G(m,m' ,5) for m,m' ^> 1 extends over a non-negligible 
region in 5 relative to the typical thermal values of m and ml . G(m, m', 5) has an integrable 
square-root singularity on its boundary in the m-m' plane. That singularity occurs when 
the point of stationary phase is the origin in position space. The stationary phase evaluation 
takes account not only of how likely it is to find the particle at the point of stationary phase 
but also how rapidly the phase begins to change as you go away from that point. It is this 
latter aspect that is optimized by applying the impulse when the particle is at the origin, 
and that is where the maxima of G(m, m', 5) occur. 
Combining these factors, one finds, e. g., for T < T c , 

j m,m' poo 
dcr Bo S e 



dtt 



poo pm 

2 / dm / dm' G(m,m',5)V(m,m',T,0) 

rp3 poo px poo 

— dx _dy zdz{(e z+x -l)-\e z+y -I)- 1 

71 Ja/A Jx+a-2y r ax JO 



a/4 J x+a— 2y/ax 

(y-x-a + 2 v ^axy 1/2 (x + a + 2^/ax - y)~ 1/2 } (10) 



where 



5 2 

a= 2T- 

Hence, da^^ /dfl is T 3 (remember that T ~ N^ 3 ) times a function / of the ratio a defined 
above. When a ~ 0(1), f(a) is 0(1). When a > 1, f(a) is 0(exp(-a/4)). And, when a < 
1, /(a) is still 0(1). To ascertain this last feature, note that the range of integrated y vanishes 
as a — > 0; y can be estimated as x in the z integral, yielding a log(l/x) integrable singularity 
upon integrating z, while the singularities in the y integral itself are also integrable, yielding 
approximately 7r/2. Hence, the apparent divergence of the z integral (as x, y — > 0) is not 
realized, even for a — > 0. A summary of these estimates can be represented as 

dch. 



~ N P e~sr 



Bose 

dn 

for all 5. 

This is the crudest of the explicit determinations of the coherent processes (simply be- 
cause eq. ( fLC|) is the furthest from a closed form in terms of elementary functions). However, 
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it is also the only coherent process whose differential rate never exceeds the Rayleigh rate by 
a divergent factor as iV — > oo. So there is no range of angle or temperature for which it is the 
overwhelmingly dominant process. A numerical determination of f(a) is certainly feasible, 
but the detailed result would shed no light on the present, general discussion. It would be 
needed, though, in the context of a particular trap potential for a particular experiment if 
one wanted to account for all the light at each angle to better than a factor of two or three. 

From the estimated behaviors given above, one can estimate the m, m! contribution to 
the cross-section: 



which is down by a factor of log(T) relative to o" B '™ e , at least when N is comparable to N e , 
because the differential m <-> m! cross-section does not grow like T/5 2 as 5 2 decreases below 
T. 



I have chosen units and notation to emphasize the angle, temperature, and number 
dependence of light scattering off an optically thin cold atom trap. The angle is best repre- 
sented in a general discussion by the momentum transfer, S. I have focussed on temperatures, 
T, near the Bose-Einstein condensation temperature, T c . And I have used semi-classical an- 
alytic techniques to evaluate the leading large- N (particle number) behavior of each type of 
contributing process. 

Rayleigh scattering or ordinary fluorescence corresponds to the • 1 term in the factor 
riiirif + 1) that characterizes boson transition rates. Hence, the sum over all possible initial 
and final atom trap states gives a factor of N for the Rayleigh rate. Because all final states 
are weighted equally, there is no dependence on 6. 

Diffraction reflects the particle density. The characteristic 5 of each sub-process is the 
inverse of the density's characteristic length scale. In a harmonic potential, using natural 




V. SUMMARY AND CONCLUSIONS 
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units, the ground state has size 0(1). Thermally occupied states have size C(T -1 / 2 ). Hence, 
the most prominent (largest angle) feature of diffraction occurs for T <T C due to scattering 
off the condensate. The condensate diffraction rate is proportional to Nq for the region 
S < 1 (in the natural oscillator units used throughout). For larger 5, this term drops off 
with an exp(— S 2 /2) because there just is no more probability to find condensate particles 
at larger momenta. The relevant rif is, indeed, rii — 1. One must also, in principle, add 
all possible struck atoms coherently. However, the diffraction from the diffuse, thermal 
cloud corresponding to excited trap states has support only for 5 — > as N — > oo. Because 
diffraction is limited in S, in contrast to the (^-independent Rayleigh scattering, the diffraction 
contribution to the total cross section has a factor of (5/ki) 2 , where 5 is the mean momentum 
transfer for the process and ki is the initial photon momentum. 

"Bose-stimulated scattering" is used here to refer to the rii • rif term, where coherence 
among the rif bosons in the final state leads to an enhancement of the rate relative to the 
totally incoherent Rayleigh rate. This is a dramatic feature of gases of Bose condensed 
atoms for T < T c when several of the lowest lying states have large occupation numbers. 
The characteristic 5 is of order the typical thermal particle momentum. This momentum 
is 0(T 1 / 2 ). However, the sum over i and / does not give cross-sections proportional to 
N N e or N 2 , where N e = N — N , even when integrated over all possible <5. Instead, it 
was shown explicitly that only 0(N^ 3 ) of the N e particles in excited states have an overlap 
with the ground state after absorption of the photon recoil momentum. The reason that the 
rii (or rif) factor does not lead to an N e in the final result is that the typical excited state 
particle has an orbit or wave function that is inherently three dimensional. The 1/3 power 
comes just from the reduction from the three-dimensional ensemble to the one dimension of 
contributing, active participants. Transferring a particle from one typical occupied excited 
state to another with a single impulse from the photon proves to be, if anything, a bit harder. 
Again, only 0(N^ 3 ) of the N e possible final particles contribute. 

The scattering rate from the <-> m processes grows like T/5 2 for angles corresponding 
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to 5 2 < 2T, and this T/5 2 is the enhancement factor over the Rayleigh differential rate in the 
corresponding angular region. The angular size of this region grows with N, since T m N^ 3 . 
This is in contrast to diffraction from the condensate, which is limited to 5 ~ 0{1). The 
m <-> w! rate, though monotonically increasing with decreasing 5, is never substantially 
bigger than the Rayleigh rate. Also, both Bose-stimulated processes fall off exponentially 
with 5 2 for 5 2 > 2T again because of the absence of particles with higher momenta. 

The heart of the calculation of the Bose-stimulated processes is the estimate of the 
efficacy of a photon momentum transfer, <5, to effect a transition between occupied trap 
states. The details of the transition matrix elements will depend on the explicit trapping 
potential. Whatever the potential, the integration over all S is two- rather than three- 
dimensional because of the constraint of overall energy conservation. The implies that even 
integrating over all kinematically accessible 6 cannot produce a Bose stimulation factor that 
is a significant fraction of N. Instead, it is generically a fractional power. 

If one observes the photon energies with high resolution (instead of integrating over 
them as done in the calculations presented here), there will be a strong angular dependence 
to the line shape. In particular, the different processes discussed here produce different 
photon spectra. At the smallest angles, the dominant diffraction leaves the photon energy 
unchanged. The Bose-enhanced scattering, dominant at somewhat larger angles, occurs with 
photon energy shifts of 0(T), while Rayleigh scattering has a characteristic, larger shift due 
to atom recoil. 
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